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We have used first-principles calculations to study the structural, electronic, and thermodynamic 
properties of the three known forms of Rh-S chalcogenides: Rh2S3, Rh3S4, and RhiySis. Only the 
first of these materials of interest for catalysis had been studied previously within this approach. 
We find that RhirSis crystallizes in a Pm3m centrosymmetric structure, as believed experimentally 
but never confirmed. We show band structures and densities of states for these compounds. Finally, 
we also investigated the ground state structure of solid sulfur (a-S), one of the few elements that 
represents a challenge for full first-principles calculations due to its demanding 128-atom unit cell 
and dispersion interactionis between units. 

PACS numbers: 61.50.-f 71.20.-b 
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There is currently a growing interest in rhodium sul- 
phide materials because of the potential of transition- 
metal chalcogenides as catalysts in fuel cells, substituting 
expensive and scarce platinumj-'- Rhodium sulfide pos- 
sesses promising characteristics, from significant activity 
in acid electrolytes to selectivity towards oxygen reduc- 
tion reaction in large concentrations of methanol.'^ Re- 
cently, a commercially available carbon-supported Rh^jSj, 
electrocatalyst has shown a superior stability and perfor- 
mance over a sample of other chalcogenides, while being 
free of the toxicity problems intrinsic to those contain- 
ing selenium.^ Other than in fuel cells, rhodium sulfide is 
useful as a catalyst in situations like efficient electrolysis 
of aqueous HCl solutions for the industrial recovery of 
high- value chlorine gasi^ 

The Rha;Sj, material used for electrocatalysis contains 
a mixture of the three phases of rhodium sulfide that are 
positively known to exist: Rh2S3, Rh3S4, and RhiySis. 
Despite of the interest on their catalytic properties, not 
much is known about these materials, in part due to their 
complex crystal structures (with 20, 42, and 64 atoms 
per unit cell, respectively). In this situation, we have 
used first-principles calculations as an effective tool to 
extract valuable information about these systems. In the 
same spirit, and as a way to validate our methodology, we 
also investigated the structural properties of the ground 
state of pure sulfur (a-S). Sulfur is one of the very few 
elements for which full first-principles relaxation calcula- 
tions have not been performed (boron being another ele- 
ment known for the difficulties that its crystalline ground 
state presents^). 

The rest of this paper is organized as follows. We de- 
scribe the method we used in Section II, present our main 
results in Section III, and summarize our work in Section 
IV. 



For our calculations, we have used density-functional 
theory (DFT)^ in the Kohn-Sham framework^, as im- 
plemented in the PWSCF program of the Quantum- 
ESPRESSO package.^ For the exchange correlation 
functional we employed the generalized-gradient approx- 
imation (GGA) functional developed by Perdew, Burke, 
and Ernzerhof (PBE)fiS since it is known that GGA gives 
somewhat better results than the simpler local-density 
approximation when describing the structural properties 
of this kind of transition-metal sulfidesJ^ The electron- 
ion interactions are described through the use of ultrasoft 
pseudopotentialsri^ These pseudopotentials were gener- 
ated with scalar-relativistic calculations in the A(f5s^ 
configuration for Rh, and in the 3s^3p'* configuration for 
S. More specific details about the pseudopotentials can 
be found in Ref. 14. 

The Kohn-Sham equations^ were solved self- 
consistently by using a plane- wave basis set, and 
the integrations in reciprocal space were performed us- 
ing Monkhorst-Pack grids The amount of plane- waves 
and grid points used depended on the system studied, 
and more details are given in the next section. 



III. RESULTS 

As a way to control the numerical approximation in- 
volved in our methodology, we studied the properties of 
the isolated crystalline elements Rh and S, and we com- 
pared the results with the well known experimental ones. 
In the process we have done what we believe are the first 
DFT calculations with full relaxation of the 128-atom 
unit cell of the most stable phase of S at ambient con- 
ditions. Then we then applied the same methodology to 
compute the structural and electronic properties of three 
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Rh-S compounds: Rh2S3, Rh3S4, and Rhi7Si5. 

A. Pure elements 

1. Rhodium 

Rhodium is a precious metal that crystaUizes in a face- 
centered cubic lattice with one atom per unit cell, corre- 
sponding to space group #225 {FmSm). 

Our DFT calculations for crystalline Rh use a plane- 
wave cutoff of 30 Ry and a 4 x 4 x 4 k-point mesh. We 
obtained a optimized lattice parameter of 3.862 A (Rh- 
Rh distance of 2.731 A) and a bulk modulus of 258 GPa. 
This shows a good agreement with experimental values 
(3.80 A and 270 GPa, respectively), and with other the- 
oretical calculations (see Ref. [H and references therein) . 

2. Sulfur 

Sulfur can crystallize in around thirty different phases 
depending on pressure and temperature, a number that 
makes it the chemical element with the highest number 
of allotropes in the periodic table. iS, At ambient pres- 
sure it forms structures with a very small packing factor. 
Twenty of those are composed of molecular rings with six 
to twenty sulfur atoms each. There also exist polymeric 
structures formed by molecular chains. As pressure in- 
creases the structures get more closely packed, and in the 
megabar regime sulfur becomes a metal with a supercon- 
ducting temperature of 17 KJi 

At ambient conditions the most stable phase of sulfur 
is a-S, in which molecules crystallize in the orthor- 
rombic space group #70 {Fddd) . Figure [T] (a) shows the 
unit cell of a-S, containing sixteen Sg molecules. Due 
to the interactions with its neighbors, the original D^d 
symmetry of each Sg molecule is lost, and instead it has 
C2 symmetry with the rotation axis along the c-axis of 
the solid, as shown in Fig.[Hb). 

We have performed DFT calculations for a-S using a 
plane-wave cutoff of 30 Ry and the F point for the re- 
ciprocal space integrals. The plane-wave cutoff was high 
enough as to provide results that were converged bet- 
ter than 1 part in 1000 for the bond length, bond angles, 
and dihedral angles of the Sg molecule. We tested conver- 
gence with respect to k-point sampling by doing a calcu- 
lation at fixed experimental lattice parameters using four 
grid points, including the F point and three more points 
at the boundary of the Brillouin zone. After optimization 
the geometrical features of the atomic arrangement did 
not show any important change, with the bond lengths 
varying less than 0.4%, and the angles varying less than 
0.8%. Table |T] shows the structural results we have ob- 
tained once the atoms were relaxed (using both in the 
fixed experimental cell vectors, and free cell vectors) and 
a comparison with experimental data. The lengths of 
the cell vectors are larger than the experimental results 




(b) 

FIG. 1: (a) Structure of a-S viewed from the [110] axis, with 
the [001] axis along the lenght of the page; the unit cell (white 
background) contains sixteen Sg molecules, (b) One of the Sg 
molecules, showing the C2 axis along the [001] direction of 
the crystal; the labels on the atoms will be used in the text 
to identify them. 



by 22.0%, 13.8%, and 6.4%. However, the S-S bonds 
in every Sg ring are better than 1% in agreement with 
experiments. Moreover, these rings show very similar 
geometrical features when the cell used is the smaller ex- 
perimental one, or when the computations are performed 
on an isolated ring. These facts are one more illustration 
of the inability of functionals like PBE to deal with non- 
local van der Waals interactions like the ones between 
Sg rings in orthorhombic sulfur. Very promising devel- 
opments within DFT to deal with systems where these 
interactions are important exist and they have been 
made computationally effective very recentlyi^ 



B. Rh-S compounds 

The rhodium sulfide compounds seem to have been in- 
vestigated first by Juza, Hiilsmann, Meisel, and Blitz, 
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TABLE I: Lattice parameters, bond lenghts, bond angles, and dihedral angles for a-S as obtained experimentally,— and 
computed in this work. The three theory columns refer to a fully relaxed cell, to a cell with fixed experimental lattice vectors, 
and to the relaxation of the isolated Sg molecule. 



Experiment Theory (this work) 

a-S a-S a-S (fixed cell) Isolated Sg 



a 




iU.4D4D A 


1 O 'TTO A 

iz.77z A 


(i(J.4D4D Aj 




b 




12.8660 A 


14.639 A 


(12.8660 A) 




c 




24.4860 A 


26.045 A 


(24.4860 A) 




S(l)-S(l') 




2.055 A 


2.070 A 


2.059 A 


2.070 A 


S(l)-S(2) 




2.054 A 


2.071 A 


2.067 A 


2.070 A 


S(2)-S(3) 




2.057 A 


2.070 A 


2.066 A 


2.070 A 


S(3)-S(4) 




2.056 A 


2.072 A 


2.069 A 


2.070 A 


S(4)-S(4') 




2.050 A 


2.068 A 


2.061 A 


2.070 A 


Mean S-S 




2.055 A 


2.070 A 


2.066 A 


2.070 A 


S(l')-S(l)-S(2) 




109.02° 


109.4° 


110.3° 


109.3° 


S(l)-S(2)-S(3) 




108.01° 


109.5° 


109.0° 


109.3° 


S(2)-S(3)-S(4) 




107.39° 


109.1° 


108.5° 


109.3° 


S(3)-S(4)-S(4') 




108.42° 


109.2° 


109.0° 


109.3° 


Mean S-S-S 




108.2° 


109.3° 


109.2° 


109.3° 


s(2')-s(r)-s(i) 


i-S(2) 


95.26° 


96.8° 


94.5° 


97.2° 


S(r)-S(l)-S(2)-S(3) 


98.12° 


97.1° 


96.4° 


97.2° 


S(l)-S(2)-S(3)-S(4) 


100.81° 


97.1° 


98.8° 


97.2° 


S(2)-S(3)-S(4)-S(4') 


98.84° 


97.2° 


98.2° 


97.2° 


S(3)-S(4)-S(4')- 


-S(3') 


96.94° 


97.6° 


96.9° 


97.2° 


Mean S-S-S-S 




98.5° 


97.2° 


97.3° 


97.2° 



who found several phases for these chalcogenides in 
1935.^^ One of them, believed to be "RhgSg" , turned out 
to be a superconductor with a transition temperature of 
5.8 K.^^ Almost thirty years later, Geller showed that 
its composition was actually RhiySis, and that this is a 
cubic material belonging to one of three possible space 
groups: Pm3,m (#221, Oh), Pi3m (#215, Ta), or P432 
(#207, O)^ It was not possible at the time to make 
a choice of the most probable of the three space groups 
to which the crystal might belong. No further attempts 
seem to have been made in order to pursue this ques- 
tion. Most of the literature published afterwards refers to 
Rhi7Si5 as the centrosymmetric structure Prruim, since 
that is the one for which lowest standard errors were ob- 
tained in the crystallographic analysis, even if Geller him- 
self writes that "this in itself does not necessarily mean 
that it is the most probable one" ^ A mineral with com- 
position Rhi7Si5 receives the name of prassoite or mias- 
site. 

Shortly after RhiySis was characterized, the struc- 
ture of Rh2S3 was also determined. This is a diamag- 
netic semiconductor with orthorhombic space group Pbcn 
(#60, D2h)- The authors also pointed out the impossi- 
bility of preparing a hypotetical "RhS2" phase that had 
been reported previously, stating that attempts to do 
so led to a two-phase mixture containing elementary S 
and Rh2S3. Rh2S3 is found as a mineral called bowieite, 
where part of the Rh atoms are substituted by Pt and Ir. 

The original work by Juza et al. described also a 
rhodium sulfide of the form Rh3S4.^^ However, it was not 



until 2000 that the characterization of this structure was 
published,^*' after Rh3S4 was synthetized at high temper- 
ature. This is a monoclinic structure with space group 
C2/m (#12, C2/1). It corresponds to a metal that shows 
a weak paramagnetism down to 4.5 K.^^ Recently, the 
mineral kingstonite of ideal composition Rh3S4 has been 
found by the Bir Bir river in Ethiopia, although again 
with Ir and Pt substituting to some degree Rh atomsj21 
In the following we apply a first-principles method to 
investigate the three materials described (Rh2S3, Rh3S4, 
and RhiySis). These are, as far as we know, all the 
rhodium sulfides that have been characterized experi- 
mentally Only Rh2S3 had been previously studied from 
first- principles 



1. Structures 

We relaxed the unit cell vectors and atomic positions 
of the three materials studied. After careful tests, we 
used a plane-wave cutoff of 40 Ry and Monkhorst-Pack 
grids of 2 X 2 X 2. This implied using eight non-equivalent 
k-points for Rh2S3 and Rh3S4, and four non-equivalent 
k-points for Rhi7Si5. 

Figure[2{a) shows the atomic arrangement we obtained 
for orthorhombic Rh2S3. The structural information for 
this material is summarized in Table |TT1 In this crys- 
tal there is only one inequivalent Rh atom, that is sur- 
rounded by six S atoms located at the vertices of a dis- 
torted octahedron. There are two inequivalent S atoms, 
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TABLE II: Lattice parameters, Wickoff atomic positions, and 
nearest-neighbor distances in Rh2S3 (space group 60, Pbcn). 
The experimental data is from Ref. |25| . 



Experiment Theory (this work) 





a 


8.462 A 


8.577 A 




b 


5.985 A 


6.049 A 




c 


6.138 A 


6.211 A 


Rh(l) 


8d 


0.10645 


0.1071 






0.2517 


0.2513 






0.0338 


0.0331 


S(l) 


8d 


0.1518 


0.1512 






0.3906 


0.3895 






0.3930 


0.3941 


S(2) 


4c 












0.9525 


0.9528 






1/4 


1/4 



both surrounded by four Rh atoms forming a distorted 
tetrahedron. The whole structure can be described as an 
alternating stacking sequence of octahedron pair layers 
along [010], where each layer is related to a neighboring 
one by a 6-glide reflection. The cell parameters computed 
agree very well with experiment, the errors being 1.4% 
(a), 1.1% (6), and 1.2% (c). Apart from slightly larger 
bonds, the PBE approximation gives excellent geometric 
predictions for this compound, finding exactly the same 
patterns of distortion that are found experimentally. The 
closest Rh-Rh distance is 3.208 A, substantially larger 
than the Rh-Rh bond in pure Rh (2.731 A). Our findings 
are very close to those of Raybaud and coworkers^i, who 
used a similar methodology. 

Rh3S4 has a complex monoclinic unit cell with 42 
atoms, of which four Rh and four S are inequivalent. This 
structure is shown in Fig. [DJb) and its geometrical prop- 
erties are listed in Table Hill together with experimental 
values. Two of the Rh atoms are still six-coordinated 
with S atoms forming a distorted octahedron. However, 
the other two inequivalent Rh atoms form Rhg clusters 
surrounded by S atoms. In these clusters, the Rh-Rh 
bond lenghts are similar to the ones in pure Rh. Again, 
our DFT results show good agreement for lattice param- 
eters (with errors of 1.9%, 2.1%, and 1.8%, for a, 6, and 
c) and reproduce very well all the features of the atomic 
distribution in the material. 

As mentioned in the beginning of this subsection, the 
crystallographic analysis was unable to conclude whether 
Rhi7Si5 belongs to group Pm3m (#221, Oh), PiSm 
(#215, Td), or P432 (#207, 0)M We have performed 
first-principles calculations starting from each of these 
64-atom unit cell structures, including also atoms ran- 
domly displaced from them. In all cases the final relaxed 
state corresponds to the cubic Pm3m case, so we suggest 
that this is indeed the ground state of RhirS'is. The 
structure and crystallographic data obtained are con- 
tained in Fig. ^c) and Table IIVI The four inequivalent 
Rh atoms have different enviroments; the central one at 




(c) 

FIG. 2: Atomic structures obtained in this work for Rh-S ma- 
terials: (a) 3x2x3 supercell for orthorhombic Rh2S3, showing 
one of the two layers of face-sharing RhSe octahedra present 
in the primitive 20-atom unit cell; (b) 2x2x2 supercell for 
monoclinic Rh3S4, showing the octahedra surrounding two of 
the four inequivalent Rh atoms, and the bonds between the 
other two types of Rh atoms; (c) 2 x 2 x 1 supercell for cubic 
RhiySis, showing the octahedra surrounding one of the four 
inequivalent Rh atoms. 
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TABLE III: Lattice parameters, WickofF atomic positions, 
and nearest-neighbor distances for Rh3S4 (space group 12, 
C2/m). The experimental data is from Ref. '2fj, 



Experiment Theory (this work) 





a 


10.29 A 


10.491 A 




h 


10.67 A 


10.899 A 




c 


6.212 A 


6.326 A 






107.7° 


107.9° 


Rh(l) 


»j 


0.36278 


0.3616 






0.14491 


0.1464 






0.4501 


0.4499 


Rh(2) 


4i 


0.3509 


0.3495 
















0.0546 


0.0550 


Rh(3) 


2c 






















1/2 


1/2 


Rh(4) 


45 












0.1596 


0.1611 












S(l) 


4i 


0.4165 


0.4157 
















0.7291 


0.7256 


S(2) 


8j 


0.1262 


0.1260 






0.1582 


0.1586 






0.3854 


0.3853 


S(3) 


4i 


0.1182 


0.1171 
















0.8883 


0.8883 


S(4) 


8j 


0.3493 


0.3501 






0.2111 


0.2089 






0.0990 


0.1004 



(1/2,1/2,1/2) is surrounded by six S atoms forming a 
perfect octahedron, while the other three Rh atoms have 
four S neighbours each. The Rh-Rh bonds get stronger 
than in the previous two materials, being around 3.5% 
shorter than in the case of pure Rh. Our DFT calcu- 
lations reproduce all the experimental structural details, 
showing a lattice parameter and bond lengths around 1% 
larger. 



2. Electronic properties 

The band structures and densities of states computed 
using first-principles calculations can be seen in Figs. [3] 
and [31 As the Rh content increases, the materials go 
from narrow-gap insulator (Rh2S3) to metals (Rh3S4 and 
RhirSis). 

The band gap we obtained for Rh2S3 is 0.10 eV, in 
excellent agreement with the results of Raybaud and 
coworkersi^, who obtained a gap of 0.09 eV. Our density 
of states is also very similar to theirs; the Ad electrons of 
Rh are split in the t2g and Cg manifolds by the octahedral 
field, with the Fermi level laying in between. 

The higher concentration of Rh in Rh3S4 pushes some 



TABLE IV: Lattice parameters, Wickoff atomic positions, 
and nearest-neighbor distances for RKitSis (space group 221, 
Pm3m). The experimental data is from Ref. |24|. 



Experiment Theory (this work) 



a 9.911 A 10.025T 

Rh(l) 24m 0.3564 0.3574 

0.3564 0.3574 

0.1435 0.1426 
Rh(2) 6e 


0.2338 0.2373 
Rh(3) 3d 

1/2 1/2 
Rh(4) lb 1/2 1/2 
1/2 1/2 
1/2 1/2 

S(l) 12j 0.1696 0.1696 

0.1696 0.1696 
1/2 1/2 

S(2) 12i 0.2310 0.2293 

0.2310 0.2293 

S(3) 6/ 1/2 1/2 
1/2 1/2 

0.2643 0.2630 



of the Rh atoms to form bonds of similar length to the 
ones in the pure metal. This causes the gap to close, 
being populated mainly by the d states of the Rh clusters. 
Our results show that the electronic analysis made by 
Beck and Hilbert^^ using extended Hiickel calculations is 
quantitatively correct. 

In the case of RhiySis, only one of the four inequiva- 
lent atoms is surrounded by S atoms forming a (regular) 
octahedron. There are strong bonds between the rest of 
Rh atoms, and the band structure and density of states 
for this material shows its clear metallic character. 



3. Formation enthalpy 

We have computed the cohesive energy of the three 
Rh-S materials by comparing their energies in the ground 
state with the energy of the isolated atoms of Rh and S. 
Our results are: 4.73 eV/atom (Rh2S3), 4.83 eV/atom 
(Rh3S4), and 5.05 eV/atom (RhnSis). The first of 
these is in good agreement with the experimental (4.61 
eV/atom) and theoretical (4.84 eV/atom) values re- 
ported in Ref. [111. 

A useful reference for alloys is their enthalpy of forma- 
tion, computed as the difference between their cohesive 
energy and the (weighted) one of the pure elements that 
form them. Figure contains the formation enthalpy 
for each of the rhodium sulfides, using a-S and FCC Rh 
as the end components. It shows a convex hull, as ex- 
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(c) 

FIG. 3: Band structure for (a) Rh2S3, (b) Rh3S4, and (c) 
RhirSis. The Fermi level is at eV. 

pected since the three materials are stable. In fact, as 
mentioned in the introduction, the Rh-S electrocatalyst 
commercially available contains a mixture of the threei^ 



IV. SUMMARY 

In this article we present a first-principles study of the 
three rhodium sulfides that are positively known to ex- 
ist. We have obtained the relaxed structures for these 
materials of interest for catalysis. For Rh2S3 and Rh3S4 
the crystallographic structures we find are very similar 




Energy (eV) 

(c) 

FIG. 4: Density of states for (a) Rh2S3, (b) Rh3S4, and (c) 
RhirSis. Black lines represent the total density of states, and 
orange (yellow) shadow areas represent the density of states 
projected on the d orbitals of the Rh atoms (the Rh atoms 
bonded to other Rh atoms). The Fermi level is at eV. 




0.2 0.4 0.6 0.8 
Composition (fraction of S atoms) 



FIG. 5: Formation enthalpy for Rh-S materials. 



to the ones found experimentally, while for Rhi7Si5 we 
identify one of the three configurations proposed exper- 
imentally as the atoms adopt. We have computed the 
band-structure and density of states for the three com- 
pounds, something that had only been done before for 
Rh2S3 using first-principles. We present too an analy- 
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sis of the energetics of these chalcogenides, showing the 
convex huU graph for these alloys. As part of this study, 
we also performed what we believe is the first full op- 
timization of the 128-atom unit cell of a-S using DFT, 
showing how the PBE approximation gets correctly the 
S-S bonds in the compound, but that it is necessary to 
include a further correction to take into accound the van 
der Waals interactions between 5*8 rings in the structure. 
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